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Abstract. We continue our investigation [l] of a model introduced by Adkins 
and Nappi [2], in which ui mesons stabilize chiral solitons. The aim of this 
article is to show that the energy associated to cquivariant solutions does not 
concentrate. 



1. Introduction 

Physicists worry a great deal about the divergences of quantum held theories: 
realistic theories make accurate predictions only after a mysterious process of renor- 
malization that removes divergences. Much of the motivation for new theories of 
elementary particles and gravity is the search for more fundamental theories which 
are finite. 

It is perhaps less appreciated among physicists that classical physics often leads 
to divergences (or singularities) as well. The laws of classical physics are encoded as 
partial differential equations (PDEs). They do not always have regular solutions, 
even when the boundary (or initial) conditions are regular. If a singularity arises, 
it can point to a breakdown of the theory at short distances - as it happens for 
the equations of elasticity [3j. Even if the theory does not break down, singulari- 
ties might describe new fundamental phenomena, e.g., black holes or shock waves, 
which arc unimaginable in the linear approximation. Thus, our understanding of 
a classical field theory is not complete until we know the conditions on the initial 
configuration under which regular solutions exist. 

It is not usually difficult to show the existence of weak solutions, i.e., distribution 
solutions for which the equations hold when averaged by a smooth test function. 
However, to prove that solutions are smooth ( "classical'^ is often a very challenging 
mathematical problem. Moreover, to show that an evolution equation has a local 
solution, for some finite time after the initial condition, is often much easier than 
to show global regularity. 

Even for linear equations, the understanding of the regularity of solutions came 
many decades after physicists had already moved on to other matters. For ex- 
ample, the proof of regularity for solutions of the Laplace equation needed deep 
developments in functional analysis, such as Sobolev spaces, in the 1930s. 
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*In the analysis of PDEs as well as in physics, the word "classical" has a specific meaning; but 
the meaning is completely different in each discipline. To a physicist, a "classical theory" is an 
approximation in which the quantum effects are ignored; classical equations of motion are often 
PDEs. In analysis, a "classical solution" of a PDE is one that can be represented as a function, 
not just a "weak solution" which is only a distribution. 
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The most celebrated open problem of regularity is the one for the Navier-Stokes 
equations of hydrodynamics. While this "millenium problem" has resisted many 
attempts, there has been a great deal of recent progress in understanding the reg- 
ularity of solutions for nonlinear wave equations, many of which also appear in 
physics. The simplest kind of nonlinear wave equation is 

(i) ^w + S = « 

where the metric g has signature (n, 1) and V is some polynomial. This equation 
follows from the variational principle with action 



(2) S = 



\g^d^d v ^ + V{4>) 



dx 



where the field 4> '■ H^™' 1 — > K™ 1 takes values into an Euclidean space . 

Regular solutions for smooth initial data exist |H [5] when V is bounded below 
and its degree is not too large, 

In 

(3) degV < 



The inequality is saturated when the energy is scale invariant. Thus, V(4>) ~ 4> A is 
regular for spaces of dimension n < 4. When the dimension is larger, the theory 
is said to be supercritical and the solution is expected to develop singularities even 
for smooth initial data. For dimensions which are less than the critical one, the 
proof of regularity is substantially simpler. 

It is interesting to note that a <fi 4 scalar quantum field theory is perturbatively 
renormalizable precisely when n = 3 (the spacetime is four dimensional); i.e., one 
dimension lower than the condition for the classical theory to be regular. This is 
the condition for the action (not the energy) to be scale-invariant. Again, when 
the dimension is smaller than the critical one, the proof that the divergences of the 
Feynman diagrams can be removed is substantially simpler: there are no primitively 
divergent diagrams and the theory is super-renormalizable. It is only in this case 
that we have a mathematically rigorous construction of quantum field theories [6]. 
Thus, there appears to be some deep connection between the renormalizability of 
quantum field theories of some given space-time dimension and the regularity of 
the solutions for its classical equations of motion in one dimension higher. 

An example of a renormalizable theory is the linear a model [7] with four com- 
ponents (a, it i , 7T2 , its) and 

(4) V(a,7r)= ^(a 2 + tt 2 + it 2 +tt 2 - a 2 ) 2 + 6<t 

The fields a and it represent certain particles called mesons. The a meson is 
much heavier than the tt mesons, so that an approximation in which the pions are 
massless and the a is infinitely heavy is reasonable. As it happens, this is the limit 
of e — > and A — > oo, keeping a fixed. It turns out that this limiting process is in 
fact a better description of nature than the linear a model. We observe tt mesons 
everyday, whereas a mesons, if they exist at all, are very short-lived resonances. 

In this limit we get the "nonlinear a model" for which the field is constrained to 
lie on a three-sphere of radius a. If we introduce some coordinates on the sphere, 
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the limiting equations can be written as 




g^d^d^h^dg 



where hij is the metric of the sphere and dg = \J—g dx is the usual invariant measure 
in Minkowski spacetime. The Euler-Lagrange equation of this action is called the 
"wave map equation" in the mathematics literature: it is a generalization of the 
geodesic equation to several indepedent variables, as well as the analogue of the 
harmonic map equation for a domain with Lorentzian signature. 

The nonlinear a model is not renormalizable in four space-time dimensions (the 
critical dimension of its spacetime is two). Interestingly enough, regularity for the 
solutions of the wave map equation was shown in 2+1 dimensions 8, 9, 10 . This 
is in keeping with our earlier observation that the critical dimension of a quantum 
theory is one lower than that of the classical theory. In particular, we should expect 
that the wave map equation in four space-time dimensions has singular solutions 
even for regular initial data. This has been believed by physicists for a long time, 
based on scaling arguments. Specifically, the initial configurations that contain a 
"topological soliton" evolve to a singularity in finite time. Since our paper concerns 
a possible solution to this pathology, let us digress to explain this point further. 

To have finite energy, any initial data for the wave map equation must tend 
to a constant at spatial infinity. Identifying the points at infinity, this defines 
a continuous map from § 3 to S 3 . An initial configuration of topological degree 
not equal to zero cannot be continuously deformed to a constant: it represents 
a new kind of excitation for the theory, fundamentally different from 7r mesons, 
which are just infinitesimal perturbations. For small times, such a topologically 
nontrivial configuration will move like a particle of finite size. It is of much interest 
to consruct a theory in which such topological solitons exist and are stable (i.e., 
free of singularities). Alas, the nonlinear a model (the wave map equation) is not 
that theory: scaling arguments show that it is possible to reduce the energy of a 
topologically nontrivial configuration by shrinking it. As the size goes to zero, the 
energy goes to zero. Thus, an initial configuration of degree one should evolve in 
finite time to a singularity, where the degree changes discontinuously to zero. This 
is by now well established, indeed an explicit solution with such a singularity can 
be exhibited [TT] : 



But wc should not give up on the wave map equation in three space dimensions 
entirely. The same physics that predicts the collapse of the soliton also suggests a 
way out: change the theory so that there is a short-distance repulsion which will 
halt the shrinking of the soliton at some finite size. Skyrme [12] proposed a model 
that has exactly this feature. He further made the revolutionary suggesstion that 
the topological solitons of the nonlinear a model are atomic nuclei: the degree of the 
map is precisely the total number of protons and neutrons, a conserved quantity. 
These strange ideas have been found to be valid as an effective model [131 H] of the 
fundamental theory of nuclear interactions, Quantum Chromo-Dynamics (QCD). 
The profound connections between the Skyrme model and the topological anomalies 
of QCD, discovered by Witten [T3], were crucial in establishing this relation. An 
exact equivalence to QCD can be proved in two space-time dimensions [T5] , 



(6) 
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Physics suggests that the wave map equation modified by Skyrme's repulsion will 
have regular solutions even for initial data that are of non-zero degree. However, 
proving this is a difficult problem: the Skyrme equation is even more nonlinear 
than the wave map equation. More precisely, it is no longer semilinear, but is 
quasilinear. This means that the coefficients of the top order terms are functions of 
the field, and, as a result, not even the local existence theory follows from standard 
theorems, but has to be proven directly. 

There exists however a simpler alternative to the Skyrme repulsion, suggested by 
Adkins and Nappi [5] , which leads to a semilinear equation. Physically, it describes 
the interactions of three 7r mesons with an uj vector meson field, which provides 
the repulsive force. Scaling arguments, as well as the existence of a static solution, 
suggest that this equation has regular solutions even for large initial data. We thus 
conjecture that the Adkins-Nappi model is no longer supercritical, i.e., solutions 
starting from finite energy smooth initial data are global in tim^j 

By now there is a well-established approach to proving such regularity theorems, 
which has worked for the wave map equation in lower dimensions. Due to the fact 
that our equation is semilinear, we know right away that for any smooth data there 
exists at least a local solution. First, one proves that energy does not concentrate 
(i.e., the amount of energy within a ball tends to zero as its radius goes to zero) 
and that the solution remain continuous at the first possible singularity. If we can 
also show that small energy implies regularity, we would have achieved our goal of 
proving regularity. 

Earlier [T] we succeded in showing that a spherically symmetric (equivariant) so- 
lution is continuous. In this paper, we complete the first part of the above scheme 
by showing nonconcentration of the energy, again, in the case of a spherically sym- 
metric solution. As time evolves, some energy is radiated away and what remains 
should settle into a static solution. 



2. Main results 
The model we are studying is given by the action 

(7) S = \j 9^ d„P hiM) dg+\J F^F^ dg ~ J A^f dx 
where: 

• <j) : (R. 3+ ,g) — > (S 3 ,/i) is a map from the Minkowski spacetime, with g = 
diag(— 1, 1, 1, 1), into the sphere S 3 C K 4 (with the induced Riemannian 
metric h), describing the tt mesons; 

• the 1-form A — A^dx^ is a gauge potential representing the lo meson, while 
the 2-form F M „ = d^A v — d v A^ is its associated electromagnetic field; 

• j is the flux or the baryonic current, i.e., 

p = ce'"'"' dptf d a <j> k e ijk 
with e being the Levi-Civita symbol and c a normalizing constant. 



The limit of zero mass for the u and tt mesons is much simpler to study. It is expected that 
the short-distance behavior, such as singularities, are unaffected by the masses of the mesons. 
This is why we will study the massless limit of the Adkins-Nappi model. 
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The first term in ^ is the action of the nonlinear a model, while the third one 
describes the coupling of the 7r mesons with the vector meson. 

For a field <f> of finite energy, one can associate naturally its topological charge 
(winding number) 



(8) 



Q 



cfdvPdntfdrF e ijk dx 



fdx 



In the Adkins-Nappi model, the density of the gauge potential A coincides with the 
density of the topological charge Q. 

The object of our investigation are time dependent equi variant maps of winding 
number 1 associated to ([7]), i.e., 

4>(t, r, V, 9) = (u(i, r),fl>, 9) A(t, r, ^ 0) = (V(t, r), 0, 0, 0) 

u(t, 0) = ll(t, 00 ) = 7T 

that yield the following PDE for u: 



sin2u (it — sin u cos it) (1 — cos 2m) 



(9) U tt - U rr - -U r 

A simple calculus shows that the energy norm 

POO 

(10) £[u](t) 



1 , n n, sin 2 u (u — sin u cos u) 2 
-{u 2 + ul) + — + y - 



2r 4 







r dr 



is a conserved quantity if u satisfies ^9|. 

Based on the fact that ([9| is a semilinear equation, for which smooth solutions 
exist at least locally, and u is a radial function, we can assume, without any loss of 
generality, that our solution starts at time T = — 1 and develops its first singularity 
at (0,0). In this context, our main result is 

Theorem 2.1. The local energy of a solution for ([9| does not concentrate: 



(11) 



lim 

T->0- 



|T| 



sin 2 u (u — sin u cos uY 



2r 4 



r 1 dr = 



This should be compared with what we were able to prove in pQ: 

Theorem 2.2. ([1 ) Ifu satisfies ([9]), then the following local energy estimate holds 

/'' T ' \( r \ o 9 sin 2 u (u — sinwcosw) 2 ! , , 
(12) lim / ( 1 - — I u 2 _ + u\ + — — + ^ '— r 2 dr = 



ITI 



where 

u + = Ut + u r U- = Ut — u r 

is the standard notation for derivatives along null directions. As a consequence, u 
is continuous at the origin with 



lim u(t.r) 

(t,r)->(0,0) 







Remark 2.3. The continuity argument uses only 

, r ' T ' (u — sinucosu) 2 

lim 



r 2 dr = 
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Remark 2.4. One notices that what we are missing in |12|) from claiming (11) 
is an estimate near the cone (i.e., the region where r w \T\) for the null direction 
derivative u_ . 



In proving Theorem 2.2 we used the conventional strategy (e.g., [5], [TB]) of 
obtaining local energy estimates by integrating differential identities over appro- 
priate past oriented domains. These were derived in turn by multiplying ([9]) with 
various conformal Killing vector fields. However, this method does not get us quite 
to where we want (i.e., proving (111), mainly because, in our case, the energy is 
not critical (scale- invariant), some of its terms being supercritical 

This is why here we adapt a technique pioneered by Grillakis 8], in which we 
work backwards in time (i.e., from positive times towards 0) and take into account 
the behavior of u, both on past and future oriented domains. We would still be 
integrating differential identities, but Grillakis's method has the extra advantage 
of working with two apparently different expressions for the local energy: 



E[u](T) = 



(13) 



o 
1 

71 



2 («? + <) + 



sm 2 u 



(a 



sm u cos u 



2 dr 



■ 2 

sm u 



2r 4 

(it — sin u cos u) 2 

2T 3 



where Kt = {(£, x)\ T < t — 2T — \x\ < 2T}. The first one is an integral over a 
spacelike ball, while the other is based on a light cone. They are equal because their 
domains form the boundary of a solid cone in spacetime. However, this equivalence 
gives better control on the behavior of energy on small balls. 



3. Notations, multipliers, and basic estimates 

We use the framework and notation of [5J , in which u is assumed to be a smooth 
solution for ([9| in the domain 

(14) D(T) = {(t,x)\ \x\ < min{2T -t,t}, < t < 2T} 
with initial data given by 

(15) u(2T-r,r) = f(r), T < r < 2T 
For < < < h < T, we denote 

D(t ,ti) = {(t,x)\ 2t Q -t< \x\ < min{2ti - t,t}} 

C(t ,h) = {(t,x)\t <t=\x\<t 1 } 

C(t ,h) = {(t,x)\\x\=t-2t <t 1 -t } 

D t0 = {{t ,x)\\x\<t } 

K to = {(t,x)\t <t = 2t -\x\<2t } 



^For a more detailed discussion of this aspect, we refer the reader to Remark 2.6 in [T] 



NONCONCENTRATION OF ENERGY FOR A SEMILINEAR SKYRME MODEL 



7 




Ctto.tj) 



Using the multiplier au t + bu r + cu, where a, 6, and c are all functions of t and 
r, we obtain that u satisfies: 
(16) 



a + b 2 
ft I u + 



a ~~ b 



u_ + a 



sin u (u — svnucosuY 



2r 4 



+ CUUt — Ct - 



a + b 



i 2 _-b 



sin 2 u (u — sin u cos u) 



at d r (br 2 ) 



+ (a t + b r ) 



sin 2 u 



2r 2 
a t +r 



2r 



+ c) u t + (—a r + b t )u t u r + 



.> ^ j i. ( ^ (it — sin u cos u)^ 



a t - b r 



□c 



CU1t r 



For (a, 6, c) = (1,0,0) we derive the energy differential identity 

._,„. n (1, o 2 \ sin 2 u (u — sinucosu) 2 \ 1 „ . , . 
(17) d t U(u 2 + u 2 r ) + —^+ ( - — 4 '-) ~ -d r (r 2 u t u r ) =0 



which, integrated over the interior of the cone bounded by K t and D t , justifies ( 13 ) 



Second, integrating (17) on the domain D(to,ti), we obtain 
(18) E(t!)-E(t ) = F(t Q ,tx) 

where F(to,ti) is the flux between the time slices D to and D tl , i.e. 



F(toM) 



V2 



C(t ,ti) 



sin u (u — sin u cos it) ^ 



2r 4 



Arguing exactly as in Section 2 of 1[, where the problem was forward in time, 
we can infer from (18 1: 
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Proposition 3.1. In the framework detailed above (i.e., u is a smooth solution for 
Q in the domain D(T)), the following estimates hold: 



(19) 
(20) 

(21) 



E(t ) < E{h), (V)0 < t < h < T 
lim F(0,T) = lim ulT.T) = 

||u||l- < C{E{T)) 



where F(0,T) = limg^o F(S, T) and C{E(T)) is a constant depending strictly on 
E(T). 



Finally, integrating (17) on the domain bounded by K to , C(t ,ti), K tl , and 
C(to,ti), and using ( 19 (-( 20 ), we deduce that 



lim E(tn) = lim — =. . 



1 



sin 2 u (u — sin u cos u) 2 

h — 

r 2 2r 4 



which allows us to reduce the proof of (| 1 1 p to the one of 
1 



(22) 



lim 



1 



sin u (it — sin it cos it) ^ 



for t\ sufficiently small. 



2r 4 



4. Local energy arguments 



We proceed by integrating (16) on the domain D(to,t\) to obtain: 



(23) 



V2 I 



G(u + , sin u, u — sin u cos u) — ^c r + c t H — u 2 



C(t M) 

H(u~ , sin u,u — sin u cos u) + ^ 



Cr — C t + - ] U 

r , 



where 

G(u+ , sin u,u — sin u cos u) 
H(u- , sin u,u — sin m cos u) 



a + b 2 



a — b 



u + + (a 



u 2 _ + (a + b) 



sin 2 u (u — sin it cos uV 

h - 

sin u (u — sin u cos w)' 



2r 4 



To go beyond the basic estimates (19)-(21) and obtain sharper ones, that would 



allow us to derive ( 22 ), the first challenge is to find multipliers, other than (a, b, c) — 
(1, 0, 0), for which the integrals on K tl , K to , and C(to, ti), are comparable with the 
energy, respectively the flux. As we control the energy, the flux, and the size of u 
through |l9|)-(21 ), the second goal would be to show that f D , to t± ^I is uniformly 

bounded for < t < t x < T. 
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We can achieve these by using (a, b, c) = (l, j, |), for which one deduces^] 



< G(u+, sinu, u — sin u cos u) < — u 



sin u (u — sin m cos m)^ 



2 ' r* 

ir)A\ 1 2 sin 2 it (u — sin it cos 
l z ^v < H(u_ , sin u,u — sin it cos it) < - itl H = h 



2r 4 



2r 4 



Cr + c t + 



< 



1 



c „ 1 

0<c r -c t + -< — 



Taking into account also (pL8|) , (|T9|), and (|21j), we are lead to 
(25) 



I 



£>(t ,ti 



< C(E(T)), < t < h < T 



Next, we want to prove that the above estimate implies energy decay. We show 
initially that one can ignore a number of terms from J D ^ tg t \ I, which are also 

uniformly bounded. For (a, b, c) = (ljfjf): 
(26) 

,2 , 



D(* ,ii) 



sin u Dc 2 u sm 2u 
(ot + 6 r )-^ -u -c— — 



< C(T,E(T)), 0<t o <h<T 



which follows from (21) and 

/ r-j-w / / +/ / -rftdr < T+ sup |ccloga;| 

JD(t ,t±) " JO J2t -r J t a Jr * 



[0,T] 



Coupling (|25J-(|26J) with 

'at _ d r (br 2 ) 
2 2r 2 



+ c I u 2 + (— a r + b t )u t u r + 



a t + r 2 d r 



a t -b r b 
1 



2 6 - a 2 



which are obtained through direct computations, and the crucial factj^] 
(27) u(u — sinucosw)(l — cos2u) >0 

we derive: 



D{t M) 



(u — sin u cos it) 2 



< C(T,E(T)), 0<t Q <h<T 



This yields the existence of a sequence (s n ) n —¥ for which 



(28) 



lim 

n— > oo 



r. 2 (it — sin u cos u) 2 



A", 



Hence, we are allowed to take to = in (23) (for (a,b,c) = (l, |, 7)) to obtain: 



(29) 



£>(0,T) JC(0,T) 



H(u- , sin u,u — sin it cos u) + ( c r — c t 









(*- 


Vc t - 








r ) J 



< T < T 



4 In what follows, the notation A < B will stand for A < CB, where C is a constant independent 
of A and B. 

5 used also in pQ 
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In deducing ( 29 ) we relied also on the trivial observation 



(30) 



lim 



c(o.t) r 



, = lim 
t->o 



for / e L c 



Using the decay of the flux ( 20 1 for the integral on C(0, T) and a slight variation 

sin 2 u Dc n u sin 2u 



of (|26b, i.e., 



lim 

T-j-0 



D(0,T) r 



= 



we can infer from ( 29 1 that 



Proposition 4.1. Under the hypotheses of Proposition 3.1 u satisfies the energy 
decay 



(31) 



lim 



1 



2 sin 2 u (u — sin w cos w)^ 







Remark 4.2. One needs to compare (31) iot£/i (12) (obtained in [TJj and realize 
that (31) is a?so enough to clinch the continuity argument. 



We use now (23) for (a,b,c) = (1, 1, to derive 



■V2 



(u — sin u cos u) 2 wsin2u u(u — sin u cos it) (1 — cos2it) 



D(t ,ti) 



2sin 2 u (u — sin u cos u) 2 

h 

1 A 



C(t ,ti) 



Relying on (27), (31), and that usin2it > near the origin, one deduces 
. . f (u — sin u cos u) 2 | «sin2it | u(u — sinucosu)(l — cos2u) 



< +oo 



>D{0,T) 1 ' ' 

for T sufficiently small. 

We can then integrate (16), for (a,b,c) — (1,1 i), on the domain D(to,ti), 
bounded by C(to,ti), K tl , and C(0,ii) 



2*1 
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to obtain 

2 sin 2 u (u — sin it cos w)^ 



/ -/■ <-I 

JC(0,ti) JC(to,ti) JK 



2 



'C(O.ti) JC(t M) JK{t ,t x ) r 

r- f (u — sinucosu) 2 usm2it u(u — sinwcosu)(l — cos2u) 

+ ^2 / 5 L + — 3— + - 

JD{t ,t{) r r r 

where 

K(t ,h) = {(t,x)\h<t = 2t 1 -\x\<t 1 +t } C K{h) 



The integrals on i^(to,ii) and D(to,t\) are handled through ([19| and (32), and so 
we can imply 

lim / u 2 + = I u 2 + 

for t\ sufficiently small. 

A similar argument for (a, 6, c) = (1,5, gr) leads to 

3 n sin 2 u (u — sin u cos u) 2 
lim / -u 2 , + ' y 



^oJ d{toM) A + 2r2 4r4 



3 2 sin 2 u (u — sin m cos m) 2 



Combining these last two limits, we derive (22 1 and thus finish the argument. 
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